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r NONCONFORMING FINITE ELEMENT METHOD FOR THE SOLUTION OF 
RADIATION TRANSPORT PROBLEMS 

K All SCOK KANO* 

Attract. The straulati wi of radiation transport in the optically thick flux limited diff'iatai regime ha# 
hern Identified as one of the most time-consuming tasks within Urge simulation rode* Due to mult imat trial 
complex geometry. the radiation transport system must often he solved cn unstructured grids. In this paper, 
wt investigate the behavior and the benefits o i the unstructured T noaronforroing finite element method 
which ha ntoven to be flexible an*i effective on related transport problems, In solving unsteady implicit 
nonlinear radiation diffusion problem using Neat on and Picard linearization method*. 

Key words, nonconform^.* finite elements, radiation transport, inexact Newton linearization multigrid 
preconditioning 

Subjarf rinse Meat ion Applied and Numer cal Mathemaf v * 

I. Introduction Radiation t ransporr in aftrophynraJ phonon en a and inertial confinement fusion is 
ofun modde-d using a diffusion approximation [12 IT. 16. 2fi. 21. 22. 24] When the radiation field U not In 
thermodvruumc equilibrium with the material a coupled set of rime dependent diffusion equations is used to 
describe energy transport These equation# are highly nonlinear and exhibit multiple time and space Kales 
Implicit integration methods are desired to overcome time step restrictions 

Moor ..nmng finite-element methods have proven flexible and effective on incompressible fluid flow 
problem* turb as mcrar.prrssibh' Stake* and Navies Slokes equation# [10, II] In the P\ nonconforming 
method, the degree# of frwdorr. lie on midpoint# of edges Therefore the number of connection# of degrees 
of freedom with each others at nv*t four four at interior edges and two at boundary edges; which » the 
same vjinbu of cormerfions of degrees of freedom in structured fimt? difference methods In contrast 
Is the Pj conforming method riie number of comertiom of degrees of freedom is a* Vumi fair exrept at 
bousdarv point# and depend# the triangulate ar»d position of point* The number rf connections of 
degrees td freedon. de*ermir ^0 ihe n umh er cf nonzero entries of generated mitrices and play# an essential 
roie hi performance of para IH knpleme mat ion# because of the communication required tn kernel operation# 
like metric- vector multiplication P% nonconforming method* general, matrices tha» have a constant irnall 
nunr.bt* of rmaaro entries f t*r each rem end tbe.>fore have K<m» advantage# in parallel implements**! and 
performance. 

Because many nonlinear slliptir problem# are well solved b> conforming fio.tr rlesomt method* non- 
conforming methods are still rar* in suet) problems Ho»e\«v nonconform trig methods may resolve feature# 
of solution 4 of nonliticar pr.i blrn s not w* it rrper*»nf/d In cmlrming method# In this research, a nonron 
forming method*# •bown tn resolve v#t> sharp change# of energies on heterogeneous d.*nain# Tl*c iriult# 
are very stonier to the solutions of the finite v.4ume method with an edge-biwH flu* lim.w 19] 

To solve nonlinear problem#, one usually envoys lineanzatkn ubniquc* Vlanv modeler# use Pxvd 
®iid N e wton wthodi to bre-ar #e Picard V rcv*h'*J .** to * y to iinderstnnd and implement, bj* converges 

•:CAfE. M^i *->f» 123C. % Ai A Lm.g^y Kwn i Ceot**. . VA .-2199 mmm II: kfcMMIraar^tu Tn. 
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dtmiy HwwUm* method has a second-order convergence r at* four require* tb* Jacobian of Cite original 
•oalinnnr system In man> nseilinear problems an inexact Newton method works well, with baa storage 
«l operation count expense [f\ !r> this paper, we study the behavior of them three methods on a model 
radiation tram port problem 

beraaae the system generated from some linearization of the nonlinear problem )» usually neusymmetne 
•a aee preconditioned GMRES [23*. Aa a pracondit loner, wt consider multigrftd Multigrid represents an 
important advance In algorithmic efficiency lor the solution o' large problems [2. 3. 4. 14. 19, 25] 

To uae rault'grid. we need to define intergrid transfer operators betweeu non am forming finite-element 
apace* Due to the noo-nc»t<dn'K of rvjnmnfreming ap«m there is no natural inter grid transfer operator. In 
previous studies of the Dooconforriing multtgrid method, tlir average value nf tun adjacent elements is used to 
get the interpolated value r a nodr Nonronfont:ing multigrid with this intergrid transfer operator is a good 
solver ke linra r nTtcmr a aid -orr.e nonlinear sytvms with smooth nonlinear coefficient* [1. 5, 6 , 9, 13. 16|. 
Boirvcr thts *ntcrgrid traru'r* jperatnr do« not pre s erve positivity of f jnrtions. which is an essential part of 
radiation transport prr*»*nr.« hermit* energy and tempera* ure are always positive Therefore some nonlinear 
problem* with discontinue* » coefficients bound eorstratats on sc&tioai. and rapidly changing soluikins. 
bite the radiation transport problem, cannot use this inter grid transfer operator beODUe the cuarw level 
approximation obtained friar the (fore ievel approximation doss not sattofy solution bounds, anu one cannot 
generate the coarse level systems or solve the coarse level problems [15]. To overcome then? difficult fee. we 
nse a new and simple hrrrgr.d transfer operator that pn/'Tves positivity and solves the above mentioned 
problem However rnuhi-nd with this inrergr.d transfer operate is ilowrr than with the previous operator 
Therefore, we use the simple intergrid transfer operator to derive coarw Wd §% sterns and the average value 
inter grid operator to solve the linear system* 

The re* t of the paper it nrgan.x*d a* fallow* In tertian 2, we describe a model radiation transport 
and it* P } nonconform .ng disc* Hriat ion In ssrtion 3. we consider a disrrettzotion is time, derive the 
ImearvaMnn* by Picard and Newton method and describe the Inexact Newton method In section 4 wv 
describe pfer-*ryjir:nne d OMULS and the nonconforming nulrigrid praconditloner Numerical experiment? 
are given in section 5 . 

2. Radiation transport model and p mom .informing d Ucnt e/at ion. Under the assumption of 
an opt#r ally thick medium (short rneua free path of photons, a first pr*aapi«t statement of rad.a*.on transport 
red un* to the radial u/n diffmon hrn.t A particular idealired dir.enr m.eai form o f sveh a *yn*m. knemr. 


a* the ‘/P model can written tw 

<n.vf **. r r«. 

(22 1, 
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Here. E'x.t) reprewir* the phrton energy Ttl.tj is the material temperature r # w the r parity and 
a is the material condirtiviiy ir. tb non equilibrium <a^- tie nor. I near rourre terms «in the tight 'hand 
$k% are buoruo and govern the tnuuff r of energy bisivr, tlie riduiion field and material temperature 
Add.t*a.al rw*' linear tlir* an* generated by the parfir liar form (4 the division coefficients wh.ch are function* 
of the E and r firids In part* ui*/ me erwury diffusion coefficient. Z>,T. £> contain* the term Xt] whirl. 
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r*4an to zb* gradient of E Tnis bmrter term i§ an artificial rpftm of ensuring physically meaningful energy 
propagation speeds fie. no fitter than the speed a f light). The atomic number the material coefficient, 
wA while H may be highh variable it it only a function A portion (i.e. z * f(i.y) in two dimension*) 
The two model probkms considered in ihie atudy are taken from |19] and depicted in Figure 1. We 
monder a unit square domain of aimilar material with auanic rrjmber i * 1 and a unit square dona In of 
two dissimilar materials, where the outer region contains material with an atomic number of r * I and the 
inner region (1/3 < t < 2/3, 1/3 < y < 2/3 1 coot aim malerirJ with an atomic number of » -* 10. The tep 
and bottom walls are im listed and inlH ar.d outlet boundary* are specified using mixed (Robin) boundary 
renditions, as shown in the figure For convenience, wt represent the boundary i * 0 0 < $ < I and r • 1 . 
0 < y < 1 by Tf. and otherwise by I Ihen the boundary condition of the problem is 

E \ dE 

°" r# 

* f). on do — r ,, 

~ *0, on W 


where n m the local outward norma) vector of I he boundary 

Lquatkaw (2 2 li and (2 2 2) lorn a syren. cf coupled nonlinear partial dfflrretiual equations which 
must be discfctnco ._i and tiirt* In this w* consider a discretization in spare and will consider 

a < hscretiiauon in time in vSe next section. 

The vanetionaJ form of (2 2 1 ) and [2 2 2 can be wntter. u follows Find (£ T) £ {H : (0) n [0 T ))> 

such that 
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for all (k, r) € *// A < and for all f f 0 fmas 

We dJarrefls* 11 by using a triangular grid containing c*d*« *. ibown in Fig jo 2 The gnd is generated by 
cunnecting of thr iQid^sDta of the edges of the tr. aright flow the cn*iar«f disrret t ration T\ whkh regains 
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«df»f and 'onformt to the mafcmaJ tmorUc* brrmlmm io stick a »c% that *ntngJ# tdgee cram chia 
boundary Let h, and 7i # i 7}. lor j * 1. /, bf giem »hm J } m % parntirn of 0 info triangle* and Aj 

it Cke mapnmum diur*rt#-r of t>* elemrm* cf 7} . 

Defn.r the Pj norK'^tffxmmf finite eleminrf ppare# 

1', - {* f I J iO: rj* it Imea/ for all /C € T >f 

t is contm'jcdt ai Om midpoint? of interior edge*}. 


Then the aoneonforxing finite element diarretfranpn of (2 2 4j and (2 2 5) ran be %Tiften at . Hud 
\E*. Tk) € {Vj * (0. I#r*ajtj 2 tucti that 


j ^ v jj + ! D.(t, r»,r£. Vw4x4^ 

~ / <MT*KT* 4 - £>. »ud> - / 2yiidrt * rj, 
/n A, 

/ ^14*+ f D.'t* T, vr« r.rfy - / <T.<r»r To* = 0 . 
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for ail (a, vj € Vj and for a/! t C JO. Imaa] 

In above filiations fo pr^forw tbe integration in ppAre, ae use a three point quadrature rule <m each 
truxr> in T f Because the pfant< mHere il* drfrer* of f^dom arc defined and the quadrature point? of 
triangle are the tarn# *e ran eaeil* meripute f le int *fr* f fori on et.ch triangle and 
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for ail bam fur’/ai of l and K £ 7 ; Ala> V u i» < ron« tan* on earh frlan*> In 7 ; 

f jt a_J t € F, •* compute Vu' nr*-dH if. /), r\m :W 


3. Tima integratloD and nonlinear iteration In tills aectiuti. wv consider « dwryet iration in time 
and three nnnHaaar Iterations i * Nra-tjn fVard and inn art V-mt *> urr/*iior 



Tl* t mm demwtiwt are cherretiied aa flrtf <irder birkviH difference* »*t h lumping of the btjm matrix, 
\mAmt to an implicit arhem# which requires the solution of a nonlinear problem at aadb time vtep. This 
approach if ftrr -order accurate to time, and ii choaen merely for convenience iince the principal objective 
ia the study of the solution of the nonlinear system Higher order temporal dbcret (rations are demons* rated 
to be worth while in [IS] 

To solve the nonlinear problem (2 2.6) and (2 2.7) we consider the Pkard Uncariiation method and the 
Newton I: near trat ion method In both method*. %t nerd to solve linear syntems to get correction! at each 
noalmear iteration ate? 

The fully implicit Picard Wariravion me'hod separatee the operator! into linear parts and nonlinear 
porta and afi nonlinear parti are equated at the previous nonlinear iteration level, k - 1 Thii result! in 
the following tynem of equatm* 
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for ail ;*,e) € V). Berajae :3 3 1> and f3 3 2, are linear syrtems in '££*\7^ 4 ) f we ran easily calculate 
thesr Jacobian 

To get the oorrurtions (SE ST) in the Pica/ 1 Method at level k, we srdve the following linear syflems 
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For the M; implicit Newton line ari ration roethxS it is somewhat rr*«jr* complicated to compute the 
Jacobian at approximate solution point* 'io get the Jacobum we have to rair^let# the derivatives of the 
e>vtem with resperf u» <P, <’*) for all fu/idion- in l j * l j 
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W ^ jpV^ can U « 4#ir> f’vd'iMt *rj on enrh Wangle in 7}. 
Afr*r we have io §c4ve the linear f>tl«m« 
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fcx earh *ep »hw J* 1 » a Jarnbiari nhirh i* rnmputed by Picard’i method or Newton'* method 
In e.iher method ** need |r* robmnneM to control ihe ftep length o where 
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In tlm «twjy. we rnmrrJ (he «tep toirh by 'imply halving a until the residual of tl# updated to uimn w 
les* than the previrm* residual In thii control, w-e sr irs'iires Ui) to gr a proper step length nr w+ n<# at « 
fixed 0 tep length and perform the next noniinev Ueratton If the numb**r »d fan ires einred* a fixed number 
then w* go to next tine trpt by u»mg thr lire* approximation, whirh ha# the smallest nonlinear r^uiuaJ 

ftr.M ARK 3 1 The methoii h<u ‘uyniptatiroiln e eerond order ronttrfence for nonlmear 

Um9 and the Pteard m^tfurd ha. an/% « fax 0*4*' 'maafftnc* Homer the rttvhmg Itnear prutUnt of the 
Putrd method »» m. rt r%*\t% ooh*d U«r\ tltm* of th* S^Hn n method b^tuM the Pira+'i m*iho4 ls^$ f/u 
eontt tion Urn* u ie+sn^H ir« r*f [7J 

Tv impfnve the effirien^> rif the \Vwtuc* ltie f ryid we ran u*e ar, »»exar» Newtiai rr»»»ViH jfi Wien 
the N>wv<n tteratKin It “la: fcom convergetKe 'i> rfc^ retidua! ia la/gj i t)»er# u ru« reaaon to ard\e the 
hoear tyt^eni acrurately Hcrweyrr when »he Nen-/ir. ivriM'n i» r’ jee" i'i e «!*# reeiduai if imail' i it* 
convergence rat# of Newton * roHh >d »• (•&! tl> couplH to» »Ke nrrurarv c4 the linear eoIutK^i To e<iju»t lie 
affftount. worli done in th» l:n«wr w4>r '*& « rr#»v#-genre u>leraroe» t> employ an inexa'i >eAi<*n 


In tk •exact SrvUm approach. tW rqnuteps wc r criteria for tlx lin * *r tolwer * proportional to th* rMMtaaJ 
tt # b» n o uhae ar iteration la tgjauufi form this it 
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where >> » 10 * I0~* a the vaJjr <ued in this tend? unlex othmrbe noted We nor* that '13; ahowi how 
to adaptively trlcct yj to recover avmploticaly full record order convergence. 


4. PGMRES and multi grid preconditioning . In this section w* explain PGMRF1S, which it a 
combination cf a Krylov-baaed linear iterative method, and mult I grid which it well known at a mrrranfui 
preconditioner, at well at a acaiabfe tolvor even in unarrelerated form for n.an> pcr44crn« 

OMRES '23) It a well known toiver for non-Ftrraitian problems In prwrtice. GMRF1S ran hr restarted 
after m atepc. where m it aorr* fixed integer parameter, to tave storage bv arrr*ptmg a generally less rapid 
con verge* o* 

We describe the rettarrH PGMRLS for anhlng 


Aji • b 


(4.4.1) 


wHb preror*ditKjnmg nel/ti Bj 

PGVCRFS(nt) Algorithm 4.1. 

(I) Start CVjotr x 0 and cr> mput* r 0 * Bjib - A/x 0 >. $ * |r c '|j and rj =* r 0 /J 
(3) Itrratr For j *- 1 . , m do 

Compifi# «- K BjAjtj 
For i « I , . .) An : 

/<ii B, Vt() 

k *=* -;, ;V, 

FwMr< 

FVefipiife *f«i j * M* la **>d + 

Lot Af> 

(3) Form tie approximation toiution 
IMv «[i., >tm.. 

r ; |l ^ • / /« l 

and tH n s i» ♦ V„y„. alwff rnirwmi/rf *Vj - ff#r,yf|, y e /T 

*4) Restart 

Compute* r,* « /?yfl - Ajr m j, if latlgflrd then flop 

elar r ran put e r 0 :• 3 • tt r >« 1 and « * r */3 and go to *2) 

Aro:>ldi iteration cnurru'ti an orthogona. bajis of th" kft pr«xindtt**i*d Krylov • ubfpAce 

Span{ r > DjAjfc JBjAj)*~’T {i ) 


It user a modified Gram-Srim.idt proc*#t in width »i* n*w verier to be orthogonaJtzed !• obtained firon 
the previous \mu* in the proce>» All r*-ed«m! v»'f *r* «r»d their nr*ro# that are computed b> the algorilhn 


ovmpaod to the prarndtooned namd?. t m * Bj<l~ Ajt m ). iaaread ci the original (tmprwws- 

dkkmtd) rvftdaaJ I * A/*«. In addrtkm. cbm a no eae> acre* to the** jnpreroodRtaMd rtbdtmU, oo.Mf 
tbry art competed exphrirly So or monitor these prmfM.tinwri re* dual* to Mop PGM RES torn tan to 
nht Bmr problem 

Nort. wt consider MuhJgrid Precooditioner Hy. 

To define a multi pi d method, wv need to define intergrid transfer operator » Mw+n nonconforming 
finite siemens sparer Dr* to the non- neat edne«« of nonoofiConrJng aparev thrr* in not a natural Imergrid 
tranifer operator In previous studies of non conforming multigrid method' I 3 6 9). average vah* of two 
adjaortit element* are uaed to set the v*hjr of a node A nomonfortning mult tgrid method with this inter gnd 
trantor operator la a good solv» r f * linear svu erns and *ome nrmtt near eynerna that haw smooth nonlinear 
coefficient* 

To grt the coarse lew! apprusimatr linear «>i unr, (r* (3 3 Q). ae need coarae .evH approximation* of 
*)a nd ^ r;- 1 ) If the appr owma r e en|, it ion 7? *‘ ! ) varies rapidly *» •pare 

then some coarse level approx i r \at sons of i “ , )ma>* hew negative val net However '/£ 4 * , .TJ , * 4 “ I ) 

are retired to he prettier tar 'he computation of . Either we cannot g« negate vhr coarse lew! systems 

or they ma> become near!;. singular. making it bard to solve the cr>ane level problem# 

To overcome these diffn' .jltiea. new and simple intergrid trawler operate cuVd thr envolume 

baaed inter grid trantor operator, which preserves only piecewise constant funrt.ons [13] It « well known 
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the original as current scanning and reproduction technology 
allows. However, the following original document conditions may 
adversely affect Computer Output Microfiche (COM) rid I or print 
reproduction: 

• Pages smaller or larger then 8.6 Inches x 1 1 inches. 

• Pages with background color or light colored printing. 

• Pages with smeller than 8 point type or poor printing. 

• Pages with continuous tone material or color 
photographs. 

• Very old material printed on poor quality or deteriorating 
paper. 

If you are dissatisfied with the reproduction quality of any 
document thut we provide, particularly those not exhibiting 
any of the above conditions, please feel frae to contact our 
Directorate of User Services at (703) 767-9066/9068 or 
DSN 427-9066/9068 for refund or replacement. 
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